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Abstract. The basic ideas of second quantization and Fock space are extended to density operator 
states, used in treatments of open many-body systems. This can be done for fermions and bosons. 
While the former only requires the use of a non-orthogonal basis, the latter requires the introduction 
of a dual set of spaces. In both cases an operator algebra closely resembling the canonical one 
I is developed and used to define the dual sets of bases. We here concentrated on the bosonic case 

where the unboundedness of the operators requires the definitions of dual spaces to support the pair 
of bases. Some applications, mainly to non-equilibrium steady states, will be mentioned. 
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gnl 1. INTRODUCTION 

Second quantization has been one of the early fields of Marcos Moshinsky and his book 
>■ ' on its applications to fermionic systems [1] has formed many of us. One of us had the 

. pleasure to look with Marcos into the extensions of this formalism to non-orthogonal 

^ \ bases and the introduction of a dual basis and corresponding dual sets of operators in 

O ' this context [2]. Since these concepts have become very useful for open many-body 

systems and recently one of us has generalized them to what may be called "third 
quantization" by applying second quantization and the Fock space concept to density 
operator spaces of mixed states and spaces of observables of fermionic systems [3]. 
The use of a set of dual bases becomes essential to maintain a simple algebraic form 
for the problem. Third quantization was successfully used to explicitly and elegantly 
^ \ solve situations with quadratic (or quasi-free) Hamiltonians and linear coupling to an 

■ environment via the Lindblad operators [3, 4] or via the Redfield model [5]. The next 

step is to follow Marcos Moshinsky to the harmonic oscillator [6], i.e. to bosonic systems 
[7]. This step implies not only the introduction of a dual basis, but actually of a dual 
space. We shall give a view of these recent developments emphasizing the elegance and 
efficiency of the methods more than any specific application. We shall present basis sets 
for the dual bases in a rather formal way and then show that it is the algebraic structure 
that determines the actual implementation. The implication, which this has in terms of 
enveloping algebras, their underlying Lie algebras and the corresponding groups has not 
been explored at all. Furthermore the possibilities to use graded algebras and combining 
fermionic and bosonic degrees of freedom in such a treatment, is very enticing to any 
disciple of Marcos Moshinsky. We shall here try to lay the foundations and give a taste of 
things to come, in the hope of enticing other participants of this meeting to help explore 
this line of thinking. Proceeding in reverse order of temporal development, we shall first 



present the bosonic case, with its richer structure and then go to the fermionic problem. 
Finally an outlook will be given on what should be done to obtain a more complete 
theoretical descriptions of each of the cases and for their algebraic unification. 

2. BOSONIC SYSTEMS 

Let us consider a Hilbert-Fock space of n bosons. Elements of can be generated 
from a particular element G M' , called a vacuum pure state, and a set of n unbounded 
operators over M', ai,...,an, which annihilate v/q, aj\Vo>— 0^ and their Hermitian 
adjoints aj, . . . ,a^, satisfying canonical commutation relations (CCR) 

[ay,flj] = 5j^k, = = 0- 

Let us define a pair of vector spaces and J^', such that contains trace class 
operators, such as density matrices, and Jt' contains unbounded operators representing 

physical observables. We choose a specific space of observables J^' and define a 
subspace Jf^ of trace class operators over J^, such that p G ^ if and only if Ap is 
trace class for any A e J^'. Thus J^' and Jt' form a dual pair of Hilbert spaces and we 
will later choose a dual pair of bases, one from each of these spaces. 

For instance, we may chose J^' as a linear space of all (unbounded) operators whose 
phase space representation of the operator is an entire function on the corresponding 
2n— dimensional phase space. Then must be restricted to operators with finite support 
in the number operator basis, i.e. to operators which have a finite number of non- 
vanishing matrix elements in this basis. Such a constraint on density matrices may be too 
restrictive for certain applications. We shall show later [using eq. (9)] how this restriction 
can be relaxed by appropriately restricting J^'. 

We now proceed with an algebraic development and conveniently adopt Dirac nota- 
tion. We write an element of as ket |p) and an element of J^' as bra (A|, and define 
their contraction or scalar product to give the expectation value of A for a state p, 

(A|p) = trAp. (2) 

We use distinct types of brackets to emphasize the difference between the spaces from 
which the ket and the bra have to be chosen. 

If b is any of the operators oy, aj, then for each p e Jt and A e J^', bp,pb and Ab, bA 

are also elements of ^ and respectively. Thus we define the left multiplication maps 
b^ and the right multiplication maps b^ over by 

b^\p) = \bp), b^\p) = \pb). (3) 

The action of their adjoint on J(f' is defined by (2) and the fact that the trace is cyclic. 
Thus 

{A\b^={Abl {A\b^={bA\. (4) 
Loosely speaking, we can also say that {P^)* = b^ and (b^)* = b^. 



Next we define the set of 4n maps dyj,dy j,j = 1 , . . . , n, v = 0, 1, 

T , - L - R 

aij = aTy, aij = aj-aj. (5) 
that have the unique properties: (i) almost-canonical commutation relations 

(ii) fly ^- left-annihilate the identity operator 

(l|«vj = (7) 
and (iii) flyj- right- annihilate the vacuum pure state |po) = | VbXVbl 

«vj|po)=0. (8) 

Writing a 2r component multi-index m — {myj E Z+; v G {0, 1},7 G {1 . . .n})-^ we 
define a dual pair of Fock bases, one for each of the spaces o^^, J/^^ as 

Their bi-orthonormality (m'|m) = 5^/^ is directly guaranteed by the almost-CCR (6). 
Here and for the rest of the paper, x = {xi,X2, ■ ■ -Y designates a vector (column) of any, 
scalar-, operator- or map-valued symbols. 
The explicit construction of the bases (9) allows us to enlarge and restrict the spaces 
and J^' such as to keep duality by always restricting one when extending the other 
or vice versa appropriately. We achieve this identifying the space with the Hilbert 
space of vectors of coefficients {Om}, 3\o) — <ym\ni), and the space with the 
/2 Hilbert space of vectors of coefficients {5^}, 3 {S\— Lm'^/w(?Zi|- Then, clearly by 
Cauchy-Schwartz inequality, | tr5c7| = | Y^m^mPm\ < °° and hence Jf and J^' are dual 
in the required sense. 

The main idea of application of the third quantization is then to express the generators 
of quantum master equations, governing the dynamics of a density matrix, in terms of 
canonical operator maps fly j , fly j. For quadratic systems, for example, such generators 
are again quadratic and can be decomposed to normal (master) modes, and thus diago- 
nalized, by means of a non-unitary analogue of the Bogoliubov-de Gennes transforma- 
tion (see [3, 5, 7]). In particular, the physically interesting non-equilibrium-steady-state 
can be constructed as the right-vacuum state of our theory. 

3. FERMIONIC OPERATORS 



In [3] the entire idea of third quantization was first presented, but the maps on the 
operator spaces were not given in terms of raising and lowering operators but in terms 



of Hennitian "coordinates" and "momenta". We shall here briefly indicate, that it can be 
done equally for fermionic (anti-commuting) raising and lowering operators. Thus the 
situation is maintained as symmetric as possible to the boson case. The main difference 
will naturally be, that the two bi-orthogonal bases will span the same Hilbert space rather 
than a dual pair. 

Thus we point out an equivalent version of fermionic third quantization, which fol- 
lows exactly the steps of the present communication, but starting instead from a set of 
fermionic operators Cy,ct, obeying canonical anti-commutation relations (CAR), 

{cj,cl} = dj,k, {cj,Ck} = {c]4} = 0, (10) 

introducing the dual sets of density operators and observables, stating (2, 3,4) and 
defining the canonical adjoint fermionic maps 

cij = ct^=#, c[j = cj^-c), (11) 

satisfying almost-CAR 

and the properties (7,8). The parity superoperator ^ is uniquely defined by its action on 
the dual vacuum states, (1 1 =^ = (1 1 , ^\po) = |po), and requiring that it anti commutes 
with all the elements of the adjoint- algebra {.^,c} = {^,c^} = 0. The difference to 
the more symmetric approach [3] is that now the canonical conjugate adjoint maps are 
not the hermitian adjoint maps Cy j c\ ^ which is however of no consequence as we 
are anyway dealing with problems in which non-normal operators enter in an essential 
way. The complete bi-orthogonal bases of the density operator space J^, and its dual, 
the space of observables J^', analogue to (9), can now be labelled by means of binary 
multi indices m, myj- G {0, 1, }, 

1^) = rK^'vjr-^ipo), (m| = (iiFK^^j)"^'^' (13) 



4. CONCLUSIONS 

We have presented third quantization for bosonic and fermionic operators states, at- 
tempting a uniform presentation. The result in itself is remarkable and very useful, but 
was essentially taken form refs [3, 5, 7]. The purpose of this presentation was to show 
the structure of the formalism presented in a form which prepares the use of algebraic 
and group theoretical techniques along the lines developed by Marcos Moshinsky. Op- 
erator states were not known and non-orthogonal bases were rarely used in this context, 
but in [2] Marcos Moshinsky together with one of the authors developed the latter for 
fermionic systems in second quantization. Even there advantage was taken of the fact. 



that the algebraic structure defines the important features of the problem. Thus combin- 
ing raising operators in the dual basis with lowering operators in the original one, one 
could mimic the algebraic structure developed for standard anti commuting operators. 
This clearly carries over to the third quantized picture. For bosonic systems the situation 
is a little more involved because we deal with operators on different spaces, yet we feel 
confident, that we can develop the techniques known as dynamical algebras or spectrum 
generating algebras also in this case, though the central interest clearly is not on any 
spectrum. The entire idea relies on the point that through bases (9) and (13) we can pass 
to equivalent infinite spaces in the first case and finite ones in the second. Limiting 
algebraic operations appropriately we never see the differences or they appear in the 
occasional use of an overlap matrix or its inverse. Note that the freedom we have in the 
fermionic case is larger than in the bosonic one, as we are never in danger of leaving the 
finite dimensional Hilbert space. 

Summarizing we may say that the construction we present provides an ideal frame- 
work for algebraic or group-theoretical developments. Yet filling the framework has 
barely begun. It also seems obvious that we can construct graded algebras also known as 
superalgebras, mixing anti-commuting and commuting variables, in the present context. 
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